Neutrino Interaction with Nucleons in the Envelope of 
a Collapsing Star with a Strong Magnetic Field 

A. A. Gvozdevt I. S. Ognev f 
"t^j-". Yaroslavl State University 

o : 

O . ul. Sovetskaya 14, Yaroslavl, 150000 Russia 

<D 

00 ' Abstract 

<n : 

The interaction of neutrinos with nucleons in the envelope of a remnant of 
collapse with a strong magnetic field during the passage of the main neutrino 
flux is investigated. General expressions are derived for the reaction rates and 
for the energy-momentum transferred to the medium through the neutrino 
scattering by nucleons and in the direct URCA processes. Parameters of the 
medium in a strong magnetic field are calculated under the condition of quasi- 
equilibrium with neutrinos. Numerical estimates are given for the neutrino 
mean free paths and for the density of the force acting on the envelope along 
i-C ■ the magnetic field. It is shown that in a strong toroidal magnetic field, the 
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envelope region partially transparent to neutrinos can acquire a large angular 



acceleration on the passage time scales of the main neutrino flux. 
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^ 1 Introduction 

Collapsing systems with rapidly rotating remnants are of great interest in astrophysics. 
Such remnants can be formed during the explosions of type II supernovae [1,2], during 
the mergers of close neutron-star binaries [3], and during accretion- driven collapse [4]. 
About 10% of the released gravitational energy is known to transform into the thermal 
energy of the remnant and to be emitted in the form of neutrinos, irrespective of the origin 
of the collapse [1-5]. A solid rotating core (of a protoneutron star or a black hole) and 
a disk (envelope) rotating with an angular-velocity gradient are formed in the remnant 
of collapse on time scales of the order of a second. The compact core of typical sizes 
R c ~ 10 km, subnuclear densities, and high temperatures (T > 10 Mev) is opaque to 
neutrinos. Being less dense and hot (T ~ 3 — 6 Mev), the envelope with typical sizes 
of several tens of km is partially transparent to neutrinos [2-5] . Within several seconds 
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after its formation, the remnant of collapse effectively cools down via neutrino radiation. 
The neutrino luminosities during this period are typically L u ~ 10 52 erg s _1 [5] . (the 
passage stage of the main neutrino flux). This neutrino flux could have a strong effect 
on the dynamics of the envelope, in particular, it could trigger its ejection [6]. Detailed 
calculations show that for spherically symmetric collapse, the neutrino absorption by the 
envelope is not intense enough for its ejection [7]. However, such envelope ejection is 
possible for a millisecond remnant of collapse with a strong toroidal magnetic field (the 
magnetorotational model of a supernova explosion [1]). 

Because of its rapid rotation and large viscosity, convection, a turbulent dynamo, and 
large angular velocity gradients arise in the remnant of collapse. These processes can lead 
to the rapid (of the order of several seconds) generation of a strong poloidal magnetic 
field with a strength up to B ~ 3 • 10 15 G on a coherence length L ~ 1 km in the core of 
the remnant of collapse [8] . Having arisen in this way, the strong poloidal field can persist 
in a young neutron star for about ~ 10 3 — 10 4 years. Duncan and Thompson [8] called 
such young pulsars magnetars. There are strong grounds for suggesting that magnetars 
are observed in nature as soft gamma-ray repeaters (SGRR) [9] or as anomalous X-ray 
pulsars (AXP) [10]. An angular velocity gradient in the envelope leads to the generation 
of a secondary toroidal magnetic field through the winding of field lines of the primary 
poloidal magnetic field frozen in the rotating envelope plasma. In this case, strong poloidal 
magnetic fields of strength B ~ 10 15 — 10 17 G can emerge in the envelope of a millisecond 
remnant of collapse in typical times of ~ 1 s [1 1] . Such strong fields can significantly affect 
the envelope dynamics even if they persist for several seconds. For example, a magnetic 
field of strength B ~ 10 17 G can produce an anisotropic gamma-ray burst [12] and, as we 
already pointed out above, can trigger the ejection of the supernova envelope [1]. 

Significantly, such strong magnetic fields can arise in the envelope during the passage 
of the main neutrino flux through it (within several seconds after collapse). Therefore, it 
is of interest to investigate the possible dynamical effects that emerge when the neutrino 
flux passes through a strongly magnetized envelope. Effects of this kind are discussed 
in the literature. For example, Bisnovatyi-Kogan [13] estimated the magnetar velocity 
gained through a neutrino dynamical kick in the case where the magnetic-field strengths 
on its hemispheres differed significantly. In this case, the kick arises because the energy 
transferred from the neutrinos to the medium in /5-processes depends on the magnetic field 
strength. However, we will discuss the dynamical effects of a different nature. Because 
of the P-parity breakdown in weak processes, neutrinos are known to be emitted and 
absorbed in a magnetic field asymmetrically, which can lead to macroscopic momentum 
transfer from the neutrinos to the medium [14,15]. Thus, the envelope region filled with a 
strong toroidal magnetic field can acquire a large macroscopic angular momentum when 
an intense neutron flux passes through it [16]. 

Here, we estimate the asymmetry in momentum transfer from the neutrinos to the 
medium during the neutrino interaction with nucleons in the envelope of a remnant of 
collapse at the stage of the main neutrino radiation. This paper has the following struc- 
ture. In Section 2, we stipulate basic physical assumptions about the parameters of the 
medium, magnetic-field strengths, and the neutrino distribution function. In Section 3, 
we derive general expressions for the reaction rates and for the components of the energy- 
momentum transferred from the neutrinos to a volume element of the medium per unit 
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time in the direct URCA processes for a strong magnetic field. In Section 4, the same 
quantities are calculated for the scattering of antineutrinos of all types by nucleons of a 
magnetized medium. Numerical estimates for the parameters of the medium, neutrino 
mean free paths, and the force density along the magnetic field are presented in Section 
5. In Section 6, we compare our results with the results of similar calculations and briefly 
discuss the possible dynamical effects of the P-parity breakdown in the envelope of a col- 
lapsing star. The calculations of our basic results are detailed in Appendices A, B, and 
C. Throughout the paper, we use a system of units with c = h = k B = 1. 



2 Physical assumptions 

Here, we investigate basic neutrino-nucleon processes in the envelope of a remnant of 
collapse with a strong magnetic field at the stage of the main neutrino radiation. We 
consider the direct URCA processes 

n + z/ e =^p + e~, (1) 

p + e~ =>- n + v e , (2) 

p + z> e =^n + e + , (3) 

n + e + => p + z> e , (4) 

and the scattering of neutrinos of all types by nucleons: 

N + Vi=>N + Vi, (5) 
N + % N + (6) 
(vi = v e , v^, v T ), (N = n,p). 

Note that /3-decay is statistically suppressed under the conditions considered. A quantita- 
tive estimate of the asymmetry in momentum transfer follows from the expression for the 
energy-momentum transferred in these processes from the neutrinos to a unit of volume 
of the medium per unit time: 



dPa (i[^) = vj n^n^/d-//)^?-, (7) 



dt 

Here, drii and dnf are the numbers of initial and final states in the element of phase 
volume; fi and // are the distribution functions of the initial and final particles; \Sif\ 2 /T 
is the square of the S-matrix element for the process per unit time; and q a is the 4- 
momentum transferred to the medium in a single reaction. An important quantity that 
characterizes the process is also its rate T defined as 

r =7 /n^n^ 1 -//)^- ( § ) 

j i f 

In particular, the neutrino mean free paths can be easily determined from this quantity: 

V 
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where N v is the local neutrino number density and r*°* is the sum of the absorption and 
scattering rates for the neutrinos of a given type. 

In the reactions of neutrino interaction with matter, we separate the medium and the 
neutrino flux passing through it. By the medium, we mean free electrons, positrons, and 
nucleons with an equilibrium Dirac distribution function: 

i 

fi 



exp(Ei/T - rji) +1 : 



where rji = n%/T\ Ei and //j are the energy and chemical potential for the particles of 
a given type. We disregard the effect of envelope rotation on the distribution functions, 
because the rotation velocity is nonrelativistic even for a millisecond remnant of collapse. 
At typical (for the envelope) densities (p ~ 10 11 — 10 12 g cm -3 ) and temperatures (T ~ 
several MeV), the e + e~ -plasma is ultrarelativistic and the nucleon gas is a Boltzmann 
nonrelativistic one. 

We consider the case of a strong magnetic field, i.e., assume that the parameters of 
the medium and the magnetic field strength are related by 

m p T > 2eB > /i 2 , T 2 > m 2 , (10) 

where m v and m e are the proton and electron masses and e > is an elementary charge. 
Condition (|1U|) implies that the electrons and positrons occupy mainly the lower Landau 
level, while the protons occupy many levels. 

The quantity (JJJ) is known to be zero [17] for a thermal equilibrium of the neutrinos with 
the medium. However, we consider the envelope region where the neutrino distribution 
function deviates from the equilibrium one. In the model of spherically symmetric collapse, 
the local (anti)neutrino distribution without a magnetic field can be fitted as follows [5,18]: 

/„ = ggfcjd 

exp(u/T u -r) v ) + l 

Here, x is the cosine of the angle between the neutrino momentum and the radial di- 
rection, uj is the neutrino energy, T u is the neutrino spectral temperature, and r\ v is the 
fitting parameter. In this paper, we ignore the effect of a magnetic field on the neutrino 
distribution function. This approximation is good enough when the neutrino mean free 
path is larger than or of the order the region occupied by a strong magnetic field. In the 
model considered in [11], this region does not exceed few kilometers in size, whereas the 
mean free paths for the neutrinos of different types are estimated to be 2-5 km. Thus, 
in our subsequent calculations, we use the neutrino distribution function (JTTJ). We will 
return to a discussion of this issue in Section 5. 



3 Direct URCA processes 

We used the standard low-energy Lagrangian for the weak interaction of neutrinos with 
nucleons to calculate the S-matrix element of the direct URCA processes. Our calculations 
of the square of the S-matrix element are detailed in Appendix A. In the limit of a strong 
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magnetic field, where the electrons and positrons occupy only the ground Landau level, 
we derived the expression 



\S lf \ 2 _ G 2 F cos 2 9 c ii 3 



T 



2L y L z V 2 uje 



exp {-Q\/2eB) 



n=0 



|M 4 



n! 



2eB 



6^ + 



|M_| 2 (Q 2 

2s (n-l)\ 

71=1 V 7 



2eE 



n-l 
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|M, 



|M_ 



4(e+j?||) (l + g a ? + + 4^(u;-A;| 



4(l-^a) 



(12) 

(13) 
(14) 



Here, 5^ 3 - ) is the delta function of the energy, the momentum along the magnetic field, 
and one of its transverse components that are conserved in the reactions; is the square of 
the transferred momentum across the magnetic field in the corresponding reaction; Q\ is 
the index of summation over the proton Landau levels; e, p\\, u, and k\\ are the energy 
and the momentum component along the magnetic field for the electron and neutrino, 
respectively; TL x L y L z is the normalization 4-volume; g a is the axial constant of a charged 
nucleon current (g a ~ 1.26 in the low-energy limit); Gp is the Fermi constant; and 6 C is 
the Cabibbo angle. 

Note that we derived the above expression for 1 5^/1 2 in [16]. It is identical to that 
from [19] and [20] when the erroneous sign of the axial constant g a is corrected (for the 
formal change g a — > —g a ). However, the paper [21], where only the term n = in the 
square of the S- matrix element was calculated, contains a difference in expression ()13j) . 
To obtain the expression for S 2 from [21], the factor (1 + g a ) 2 must be substituted for 
(1 + g 2 ) in the first term of the amplitude (|T3j) . 

For convenience, we calculate the 4-momentum, dP a /dt, transferred to the medium in 
the direct URCA processes involving neutrinos (0), (J2J) and antineutrinos ©, © sepa- 
rately. Using the T invariance of the square of the S"-matrix element for these processes 
and explicit particle distribution functions for the medium, we can represent the quantity 
(J7J) as 



dP c 



dt 



d 3 k 



exp 



u/T±8r]) )f, 



— exp ( — uo/T ± 5r] 



(15) 



Here, 8rj = (/x e + fi p — fi n )/T; k a is the 4-momentum of the (anti)neutrino, and K^ v,u ' is 
the absorption coefficient in the (anti)neutrino absorption reaction defined as 



dn p dn n dn ( 



\s if \ 2 / /„(l-/ P )(l-/ e -) 



(16) 



r 1 / P (i-/n)(i-/ e+ ) 

Note that under /3-equilibrium conditions, when r] u = 5rj = (// e + \i v — f-i n )/T expression 
(|15|) for the transferred momentum is 

dP^ } r d 3 k 



dt 



(27T) 



:kn 



1 + exp [--±n„ 



(17) 
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where 5f v,u is the deviation of the neutrino distribution function from the equilibrium one. 
Thus, the 4-momentum transferred from the neutrinos to the medium is nonzero only in 
a partially transparent (to neutrinos) envelope and vanishes in a dense /3-equilibrium 
remnant core. 

We calculated the absorption coefficient (fTH|l in the strong field limit ffTTTf) by assuming 
the nucleon gas to be nondegenerate. In this case, the calculation technique is simplified 
sharply compared to the more general case. For this reason, we discuss the details of our 
calculation of the absorption coefficient in Appendix B. Here, we present only the final 
result that was previously published [16] without explaining the details of our calculations: 

^ = 9h^ eBNM ^ + 3 9 j) _ hfa _ !)) ( t + exp ( ±a _ £ 

where the terms of order ~ eB/m p T were discarded. Here, N n , N p , m n , and m p are the 
neutron and proton number densities and masses, a = (fi e — A)/T, A = m n — m p . Note 
that the derived expression is identical to that in [19]. 

For completeness, we give the expressions for the energy density dQ/dt and for the 
rates Y of the direct URCA processes (JTJ) - (jlj): 

d Q/ dt W^- ~ Nb I T*C 3 (a, T u , Vu ) J ' ^ 



r \ ,Np_ a/tJ T 3 B 2 (a)-T 3 J u D 2 (a,T v , V . 

dQ/dt j p+e -_+ n+Ue *N B e 1 -T*B 3 (a)+T*J u D 3 (a,T u , Vi 

T \ at N p i T 3 C 2 (-a,T,, Vl 



u) 



(20) 



2G 2 F cos 2 8 C ( 9 \ 



dQ/dt j p+ ^ n+e+ ^ Jv N B \ T 4 C 3 (-a,T,,^) /' ^ 

r \ T 3 B 2 (-a)-T 3 J D D 2 (-a,n, Vi> ) \ 

dQ/dt j n+e+ ^ e ^N B 6 l-T^H + mi-a.T^,) /' ^ 

where the dimensional coefficient A is defined as 

2G 2 F cos 2 r 

N B = N n + N p . 

Our introduced functions B n , C n , andD n can be expressed in terms of the integrals as 
follows: 

oo 

Z n dZ 

e Z~a _j_ I ' 

oo 

Z n dZ 



B n (a) = 
C n (a, T u , r\ v 
D n (a,T u ,r] p 



o 

oo 



Z n dZ 

(e z ~ a + l) (e ZT / T ~-i» + i j 
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The parameter 



J »=(f f^ k ) ( I ( l + ex P (^/ T - n, 



(23) 



has the meaning of the ratio of the actual local neutrino number density to the equilibrium 
one with temperature T v . 

Our calculations of the components of momentum (j!5|) transferred to the medium 
during neutrino reradiation show that the emerging radial force is much weaker than the 
gravitational force and cannot significantly affect the envelope dynamics. Thus, of interest 
is the force component acting along the magnetic field. For a toroidal field, the density 
of this force can be represented as 



1 ri z — 1 

C^urca Ja ^\T^ 



63^ + 1 
-(3(xl)-l)e- A/T B 3 (- 



N, 
N B 

1 ^-1 
2 3^ + 1 



(xl))^ + (i - (xl) v ' Q 



dt 



dt 



Here. 



(xl) 



X 2 ujf u d 3 k 



ujf^k 



-i 



(24) 



(25) 



where (xl) is the mean square of the cosine of the angle between the neutrino momentum 
and the radial direction. As follows from this expression, the asymmetry in momentum 
transfer is nonzero for two reasons: either the neutrino distribution is anisotropic {(x 2 ) 
1/3) or energy is transferred to the medium in the URCA processes (JU - © (dQ v> p/dt ^ 
0). Interestingly, for an isotropic neutrino distribution, the force density is directed along 
the field when the medium cools down and opposite to the field when it heats up through 
the URCA processes. 



4 The neutrino scattering by nucleons 



We used the nonrelativistic vacuum wave functions of nucleons with a certain spin 
component along the magnetic-field direction to calculate the 5*-matrix element in the re- 
actions of (anti)neutrino scattering by nucleons (0) and (0). Our calculations are detailed 
in Appendix C. Below, we give only the final expression for the square of the S'-matrix 
element for the neutrino scattering by nucleons: 



\S. 



if\ 



(2vr) 4 G^ (4) 



( c I + 3cI)ujuj'+(cI-cI) (kk') + 



T 2V 3 LULU> 

+2c v c a (uk\ + uj'k\\) {S + S') - 2c 2 a (uk\ - u'k\\) {S - S') + 
+ (c 2 v - c 2 a ) (lulu' + (kk')) SS' + Acl^SS' 



(26) 
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The expression for \Sif\ 2 /T for the antineutrino scattering by nucleons can be derived 
from (|2T)|) by the change k *->■ k': 

\S lf \Uk,k>) = \S if \l(k>,k). (27) 

Here, c v and c a are the vector and axial constants for a neutral nucleon current. In the 
low-energy limit [22], 

c„ = —1/2, c a ~ —0.91/2 forneutrons, 

c v = 0.07/2, c a ~ 1.09/2 forprotons. (28) 

The conservation of energy-momentum is defined as 5^ = 8^'(V + k — V — k'), where 
k— (u, k); k' = (u f , k')\ V and V are the 4-momenta of the initial and final neutrinos and 
nucleons, respectively; k\\ and k'» are the momentum components of the initial and final 
neutrinos along the magnetic field; and S and S' are the components of the polarization 
vectors for the initial and final nucleons along the magnetic field (S = ±1). 

Analysis of the kinematics of the neutrino scattering by nucleons shows that the energy 
transferred in these reactions to an element of the medium is negligible compared to the 
energy transferred in the URCA processes. This can also be verified by a direct calculation 
using the technique detailed in Appendix C. Therefore, below, we are concerned only with 
the component of the momentum transferred along the magnetic field. This component 
arises from a partial polarization of the nucleon gas in the field, because nucleons with 
different polarizations have different energies: 

w V 2 eB 
E N = m N + g N S 



2rriN °" 2mjv' 

where is the nucleon magnetic factor (g n ~ —1.91 for the neutron and g p ~ 2.79 for 
the proton). Taking into account the energy of interaction between the nucleon magnetic 
moment and the magnetic field, we obtained the following expression for the force density 
along the magnetic field during the scattering of neutrinos of one type by neutrons or 
protons (see Appendix C for the details of our calculations): 

3[r } = m^T NsNv { ( CvcM) + clT{ujl) ) ( {xl) - 1/3 ) " 

-cl(H) - 5T(u£)) (5/3 - (xl)) + 2<£j„(<a£) - 5T„<a£>) (l - <^>)j, (29) 

where 

H) = ( I ^ n fud 3 k] ( [ f u d 3 k] , (30) 



J v and (xi) are the quantities defined in the same way as in the URCA processes, and 
N v is the local neutrino number density. During the antineutrino scattering (|SjL the 
momentum transferred to the medium per unit time is given by expression (|2*9*|) with the 
formal change c 2 a — > —c 2 a \ 

^\c 2 a ) = ^r ) (-ci). (si) 
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For the Boltzmann neutrino distribution 



/„ = $„(r,x)exp (—u/T v 
expression (|2*9"j) can be simplified: 

6G 2 F g N eB 



7i rriN 



N N N u T^4c 2 a (2-(xl)) + 
+5T U /T[c v c a ((xl) - 1/3) - c 2 a {5/3 - (xl)) 



(32) 



Note also that in the envelope of the remnant of collapse, the parameters of the distribution 
functions for the neutrinos and antineutrinos of type \x and r are virtually identical [23]. 
This allows a simple expression to be written for the total (neutrinos plus antineutrinos 
of a given type) force density along the magnetic field: 



Cj("i) 



G 2 F c v c a g N eB 
71 m^T 



N N NM)({xl)-l/?) 



(33) 



As can be seen from formulas (|29~j) and (|31j ) -([33 p . there is an asymmetry in momen- 
tum transfer along the magnetic field, as in the URCA processes, either for an isotropic 
neutrino distribution ((x 2 ) 7^ 1/3) or when the neutrino spectral temperature differs from 
the temperature of the medium (T u 7^ T). Interestingly, the force densities along the 
field during the scattering by neutrons and protons are directed oppositely. This 
follows from the fact that the expression for the force density is proportional to the nu- 
cleon magnetic factor g (recall that g n ~ —1.91 for the neutron and g p ~ 2.79 for the 
proton). Under actual conditions in the envelope of the remnant of collapse, it is com- 
monly assumed that N p /Nb <C 1- We see from expressions (j2HJ) and (|3*T|) that the force 
density is proportional to the nucleon number density. Thus, the scattering by neutrons 
mainly contribute to the quantity 3hi. For completeness, below, we give expressions for 
the dominant contribution (with the terms eBjm^T and sjl '/mjv disregarded) to the 
rate (JBJ) of scattering of (anti)neutrinos of any type by nucleons: 



7i 



N N N v (ul) 



1 - J A 1 



(34) 



Recall that c v and c a are the vector and axial constants of a neutral nucleon current, 
which change only when the type of nucleon changes [see ()28|)] in the scattering reactions. 
Note that the scattering reaction rates and, hence, the neutrino mean free paths in these 
processes are virtually independent of the magnetic field. 



5 Numerical estimates 

For the force density along the magnetic field generated in the processes of neutrino in- 
teraction with nucleons (fT|)-([B)) to be numerically estimated, the parameters of the medium 
and the neutrino radiation must be specified. Recall that we consider a sufficiently dense 
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(p ~ 10 11 — 10 12 g cm -3 ) and hot (T ~ several MeV) envelope region that is partially 
transparent to neutrinos. Part of this envelope, about 1-3 km in size, is assumed to be 
filled with a strong toroidal magnetic field of typical strength B > 10 16 G. For such en- 
velope parameters and magnetic-field strengths, condition (fL0|) is satisfied and, hence, all 
our expressions are valid. Recall that we used the local nonequilibrium neutrino distribu- 
tion function fTTf while disregarding the effect of magnetic field on it. Below, we verify 
the validity of this approximation by estimating the mean free paths for the neutrinos of 
various types and by comparing them with the size of the region filled with a strong mag- 
netic field. The results from [23] were used to fit the neutrino-radiation parameters. In 
the above paper, a numerical solution was obtained for the neutrino distribution function 
at the stage of the main neutrino radiation after spherically symmetric collapse. 

Analysis of the numerical values for and (u; 2 ) from [23] shows that the distribution 
functions for the neutrinos of various types are well fitted by the parameters: 



T Ve ~ 3.3 MeV, T Pe ~ 4.5 MeV, T^ >T ~ T^ T ~ 6.6 MeV, 

7j„ e ~2.8, 7fe~2.0, 77^-77^-2.2. (35) 

For numerical estimates, we chose an envelope region with a typical density of the medium 
p ~ 5 x 10 11 g cm -3 and assumed the magnetic-field strength in it to be B ~ 4 • 10 16 G. 
The local neutrino number densities in this region are 

N Ve ~ 9 • 10 32 cm" 3 , N 0e ^ 3 • 10 32 cm" 3 , N^ T ~ N^ T ~ 2.7 • 10 32 cm" 3 . (36) 

The following values of our introduced parameters J u correspond to these number densi- 
ties: 

J, e ~0.23, J, e ~0.05, J VMiT ~ Jp^ T ~ 0.01. (37) 

Accordingly, the mean square of the cosine of the angle between the neutrino momentum 
and the radial direction is 

(Xl) - {Xl) - 0-385. (38) 

At the stage of the main neutrino radiation (within about 1-3 s after collapse), the param- 
eters of the medium vary slowly through hydrodynamic processes compared to the time 
scales of ~ 10~ 2 — 10~ 3 s on which a quasi-equilibrium is established through the dominant 
URCA processes. For this reason, the medium is assumed to be in quasi-equilibrium: 

Tn^p Fp_> n , (39) 

where T n _ p and T p ^ n are the sums of the rates of the processes with the conversion of a 
neutron into a proton and a proton into a neutron, respectively; and dQ/dt is the total 
amount of energy transferred in the neutrino processes to a unit volume of the medium 
per unit time. Recall that under the envelope conditions considered, the URCA processes 
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give a dominant contribution to the establishment of equilibrium. Together with the 
electroneutrality condition, 

N p _ eBfi e 



(4:1) 

N B 2vr 2 iV V ; 

the quasi-equilibrium equations for the medium and (j4*U)l allow only the density of the 
medium and the magnetic-field strength to be considered as free parameters. Numerically 
solving this system of equations yields 

N 

T ~ 3.8 MeV, a ~ 2.8, — p - ~ 0.07. (42) 

Nb 

For these parameters, the neutrino mean free paths @ are estimated to be 

l Ve ~ 3 km, h e ~ 5 km, l v ^ Vr ~ fe M)Cr ~ 2.5 km. (43) 

Comparison of the neutrino mean free paths with the size of a region, ~l-3 km, filled 
with a strong magnetic field shows that the magnetic field cannot significantly affect the 
neutrino distribution functions. 

Note that the neutrino scattering by nucleons gives a contribution to the v e mean free 
path comparable to the URCA processes and a dominant contribution to the v e mean 
free path. Thus, the v e , and v T mean free paths are virtually independent of the 
magnetic-field strength. 

It is of interest to compare the contributions to the momentum transfer along the mag- 
netic field from the URCA processes and from the neutrino scattering by nucleons. The 
scattering of \x and r (anti)neutrinos by neutrons [see (33)] gives a dominant contribution 
to the force density in reactions (5) and (6), whereas the reactions with electron neutrinos 
[see formula (24) under the quasi-equilibrium conditions) give a dominant contribution in 
processes (l)-(4). Therefore, the ratio of the force densities can be represented as 



-J II _ 4:C v c a g n B 5 [r] Vx ) J Vx T Vx ( T Vx 



^(urca) g 2 _ J Bg ( fl ) Y m n \ T 



(44) 



where v x = i/^, u T , and Y = N p /(N P + N n ) is the chemical composition parameter of the 
medium. By substituting numerical parameters of the medium and the neutrinos in this 
expression, we can easily verify that it is of the order of unity, although the T Ux /m n ratio 
is small. 

Under the quasi-equilibrium conditions for the medium and (jjOj) , the expression 
for the force density along the magnetic field in the URCA processes (124)1 is greatly 
simplified. Its numerical estimate for the above parameters is 

1.4 x 10- ^ ( B — ) . (45) 

cm 3 V 4.4 x 10 16 G J K ' 



Cf(urca) _ t A w t n 20 



As expected, the scattering of neutrinos of all types by neutrons gives a numerically larger 
estimate for the total [over all types of (anti) neutrinos] force density than do the URCA 
processes: 

<>('«*) „ 3 o x 10 20 dyne / B \ ( p \ 

J « " cm 3 I 4.4 x 1016 G / Ux 10" g cm" 3 J ' (4bj 
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We numerically analyzed the quasi-equilibrium conditions in the envelope region specified 
by the range of densities 2 x 10 11 < p < 10 12 g cm -3 with the neutrino parameters at a 
given density from [23]. Our analysis shows that the force density changes smoothly with 
increasing density of the medium. In the scattering processes, the force density increases 
almost linearly; in the URCA processes, the monotonic rise gives way to a decrease at 
p ~ 8 ■ 10 11 g cm~ 3 . The value of p = 5 x 10 11 g cm~ 3 that we chose for our estimates 
is actually a point with a mean force density in the above range of densities. Note that 
the force generated in the two processes is directed along the field (i.e., the effect in all 
processes of neutrino interaction with nucleons is accumulated) and is large quantitatively. 
Interestingly, under the quasi-equilibrium conditions, the force density is completely and 
mainly determined by anisotropy of the angular neutrino distribution function in the 
URCA processes and in the scattering process, respectively. 

To discuss the possible macroscopic effects of neutrino radiation on the magnetized 
medium of the envelope, below, we give an estimate of the angular acceleration that arises 
in an envelope element under the total force density of neutrino spinup: 

where R is the distance from the envelope volume to the center of the remnant of collapse. 
To obtain this estimate, we assumed that the macroscopic momentum is transferred to 
the entire envelope element. Note that this angular acceleration is large enough to spin 
up the envelope region filled with a strong magnetic field to millisecond rotation periods 
in a time of ~1 s. 

6 Discussion 

Here, we investigate the possible dynamical effects of P-parity breakdown during the 
neutrino interaction with nucleons in the envelope of a collapsing star with a strong 
magnetic field. These effects are known to disappear in an optically thick (for neutrinos) 
medium [19], where the neutrino mean free path is much smaller than the characteristic 
size of the envelope of the remnant of collapse (see Section 5). This determines the range 
of envelope densities and temperatures under consideration. 

Another important question is: At what strengths of the magnetic fields might their 
significant effect on the processes under consideration be expected? Currently, the field 
effect on individual neutrino-nucleon processes are being intensively studied. In particular, 
Leinson and Perez [21] calculated the neutrino luminosity in the direct URCA processes in 
the core of a neutron star (a strongly degenerate nucleon medium with a typical density 
P > Pnud = 2.8 x 10 14 g cm~ 3 and temperature T ~ 0.1 MeV). It is argued that the 
magnetic field significantly affects the luminosity only upon reaching an enormous strength 
B > 7 • 10 17 G. Baiko and Yakovlev [20] showed that the direct URCA processes could 
have neutrino luminosities higher than the modified URCA processes at field strength 
B > 10 16 G. Calculations of the Rosseland mean absorption coefficients for neutrinos in 
the envelope of a collapsing star indicate that fields B > 4 • 10 15 G are required to change 
these coefficients at least by 5% [24]. As a continuation of these studies, Lai and Arras [19] 
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calculated the collision integral for basic neutrino-nucleon processes in the Boltzmann 
equation for the neutrino distribution function and its first moments. Estimates of the 
P-odd terms in the moments of the collision integral show that the effects of asymmetry 
in momentum transfer along the magnetic field in the envelope of a collapsing star can 
be significant only when strengths B > 10 16 G are reached. The emergence of such (in 
order of magnitude) magnetic fields when estimating the P-odd effects is natural. Indeed, 
only the electrons and positrons at the ground Landau level contribute to the P-odd part 
of the absorption coefficients for the direct URCA processes. The significant dynamical 
effects related to the P-parity breakdown take place at a high e + e~ -plasma density at 
the ground Landau level, which is ensured by conditions 1)10)) . 

Assuming that the plasma electrons and positrons are only at the ground Landau level, 
we derived simple analytic expressions for the force density along the magnetic field ()24|) . 
rate, and energy transferred to a unit volume of the medium per unit time (j!9|) - (}2"2"|) in 
the URCA processes, as well as general expressions for the force density and the rate 
of scattering of neutrinos of any type by nucleons. The asymmetry in momentum 
transfer along the field in these processes is nonzero only in the envelope region partially 
transparent to neutrinos ((x 2 ) 1/3, T v ^ T ).. Assuming that the medium is in quasi- 
equilibrium via the dominant URCA processes specified by Eqs. (j3T?j) and (j4*U|) . we find 
the equilibrium parameters of the medium. Numerical estimates of the force density along 
the field (|45)l and (f4"E|) for the equilibrium parameters of the medium show that in the 
sum of the processes of neutrino interaction with nucleons, the asymmetry in momentum 
transfer is accumulated and this asymmetry is quantitatively large (J47j) . 

It makes sense to compare our estimates with calculations of the same quantities in 
the processes of neutrino interaction with a strongly magnetized e + e~ -plasma [25]. It is 
important to note that the force density along the field in these process is directed along 
the field and may be of the same order of magnitude as that in the neutrino-nucleon 
processes. Thus, the asymmetry in momentum transfer along the field is accumulated 
through the P-parity breakdown in all significant processes of neutrino interaction with 
the medium. Note, however, that Kuznetsov and Mikheev [25] worked with an equilibrium 
neutrino distribution function (a Fermi-Dirac distribution with a spectral temperature 
T u ). This causes the number of neutrino states to be significantly overestimated under 
actual conditions in an envelope partially transparent to neutrinos. Indeed, our introduced 
parameter J u , which has the meaning of the ratio of the actual neutrino number density 
to the equilibrium one, is much less than unity in the envelope region with the densities 
and temperatures under consideration [see estimate (|57j)]. Thus, in this envelope region, 
the force density in the processes of neutrino interaction with a magnetized e + e~ -plasma 
is low compared to its value in the processes of neutrino interaction with nucleons. 

An asymmetry in momentum transfer along a toroidal magnetic field gives rise to 
the angular acceleration ()47|) of the envelope region filled with such a strong field. This 
acceleration is large enough to spin up this envelope region to millisecond rotation periods 
in a time of ~ 1 s. This large (in magnitude) and local (coordinate-dependent) angular 
acceleration can cause a rapid change in the gradient of angular velocities in the envelope, 
which, in turn, can change the mechanism of subsequent generation of a toroidal magnetic 
field. Indeed, in the presence of an additional angular acceleration in the envelope, the 
toroidal magnetic field can grow with time much faster than a linear law [11]. This 
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toroidal-field rearrangement can affect the mechanisms of supernova envelope ejection [1] 
and the formation of an anisotropic gamma-ray burst in a "frustrated" supernova [26]; it 
can also trigger the growth of MHD instabilities [27]. However, to investigate the effect 
of neutrino spinup on the envelope dynamics and the generation mechanism of a toroidal 
field requires analyzing a complete system of MHD equations. This is a complex problem 
far outside the scope of this paper. We hope that it will find its researches. 
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Appendix A Calculating |5| 2 for the Direct URCA 

Processes 

In the low-energy limit, the effective S'-matrix element of the URCA processes CQ)-(@J) 
can be represented as [5] 

S lf = Gf ™ 6 c J d * x [*W 7a (l + g alb )^ n<jl ] [^ e7a (l + 75 )*J (A.l ) 

(we use the 7-matrices in spinor representation with a different sign of 75 [28]). Here, \l/ e , 
^v, ^na 1 -, and ^ p n J are the wave functions of the electron, neutrino, neutron, and proton 
(at the nth Landau level), respectively; 9 C is the Cabibbo angle; g a is the axial constant 
of a charged nucleon current (g a ~ 1.26 in the low-energy limit); and o and a are the 
components of the double proton and neutron spin along the magnetic field, respectively. 
We performed our calculations in a coordinate system where B = (0, 0, B) in the gauge 
= (0, 0, Bx, 0). When the e + e~ -plasma occupies only the ground Landau level, the 
wave functions of relativistic electrons can be chosen in the form 

* e - 6XP H(£t " M - E^Mll rr J^.fTB ( r , ) 
^/2eL y L, v eBJ 



U \ „ fo 



U e = ^e + m e [ „ ,[/= , (A.2 

1 -p 3 /(e + m e )U J V 1 " 



where p = (£,Pi,P2,Ps) is the electron 4-momentum and e = \fp\ + ™ 2 e is its energy. In 
what follows, 

Xn(x) = 1 j=^= H n (x), (A.3 ) 
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H n (x) are the Hermitean polynomials, V = L x L y L z is the normalization volume. 

The wave function of the nonrelativistic protons at the nth Landau level was chosen 
in the form 



= exp [-i(E p t - P 2 y - P 3 ; 
^2EpL y L z 



[ ^\v = ^B (x-^- 

eB 



Xn(v) 





,K=-i = 







Xn-l(v)J ' 



(A.4) 



where P = (P p , Pi, P2, P3) is the proton 4-momentum, E p = m p + P|/2m p + eBn/m p is 
the nonrelativistic energy of a charge particle at the nth Landau level. The wave function 
for neutrons is 



,t, exp[-i(E n t-qiX-q 2 y-q 3 z)] r—— 

= V E n + m r , 



V2EnV 







(A.5) 



where g = (E n , qi, q 2 , q 3 ) is the neutron 4-momentum, E n = m n + (q) 2 /2m n . We assume 
the neutrino to be a standard massless Dirac particle of left-handed helicity with the wave 
function 



W + W = 1, 



exp [—i(cut — k\x — k 2 y — k 3 z) 



w 



(A.6) 



where k = (u, ki, k 2 , k 3 ) is the neutrino 4-momentum. Squaring Sif, integrating over d A x, 
and summing over the particle Landau levels and polarizations yields 



\S,. 



if\ 



T 



G%cos 2 9 c n 3 , ^ ln 
^-^_exp(-Qi/2 eB )x 



X 



£ 

n=0 



I Mi 



n! 



2eP 



|M_ 



^-f (n- 1)! V 2eB 

n=l v ' N 



Ql 



n— 1 



5 (3) 



(A.7) 



Here, 5^ 3 ^ is the delta function of the energy, momentum, along the magnetic field, and 
one of its transverse components conserved in the reactions; Q\ is the square of the 
transferred momentum across the magnetic field in the corresponding reaction; and n is 
the index of summation over the proton Landau levels. The quantities |M + | 2 and |M_| 2 
are defined as 



IMJ 2 = 



1 



4m p m n 



<r'=-l 



L a(3 = Sp p u 7 a (l + 7 5 )p e 7 /3 (l + 75) 



AT< J < J 



Sp pna'1a(l + QalblPpalp ( 1 + 9alb) 



(A.8) 

(A.9) 
(A.10 ) 
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Below, we give the neutrino, electron, neutron, and proton density matrices that corre- 
spond to the wave functions: 



p v = k(l - 75)/2, pe = (P\\ + m e )II_, 
p na > = m n (l + u)U a >,p pa = m p (l + u)U a , 
k = k a j Q , u = u a j a ,p\\ = po7o - P3l3- 
n a = (l + az 7l72 )/2 



(A.ll 



Here, IT^ is the projection operator of the density matrix for a charged fermion in a 
magnetic field (recall that a, a' = ±1 are the components of the double proton 
and neutron spin along the magnetic field). To ensure the covariance of our calculations, 
we inserted = (1, 0, 0, v), the 4- velocity of the medium along the field, in the nucleon 
density matrix (at the end of our calculations, we assume that v — 0). 

In the subsequent calculations, we used the following properties of the projection 
operator: 



u 2 a = n CT , n CT n_ 



o, n CT 7 Q || = 7 Q |in (T , n CT 7 a _L = 7 Q _ L n_ (T , = 7 5 n CT , (A. 12 



where 7 a ii = 7i j2 and •y a ± = 70,3. Since the velocity vector v of the medium has only 
longitudinal components, the nucleon spur (|A.1Q |) can be reduced to 



ly af3 



m p m r . 



i n Sp [n CT /7 a n CT 7 /3 (l - g 2 a 
+ m p m n Sp [n (7 /7 a n (7 W||7 ( 8U|| (1 +g 2 a - 2g a j 5 



(A.13 



In the subsequent calculation of this expression, it is convenient to separate out the 
contributions from identical and different nucleon polarizations. The matrices 7^ and 
■jp have only longitudinal components (j a = j a ±, 7/3 = J/3±) when a = a' and only 
transverse components {pf a = 7 a j_, 79 = r Yp±) when o = —a'. This is easy to obtain from 
the properties of the polarization operator (jA.12 |) . To calculate the nucleon spur requires 
the expressions 



Sp 
Sp 
Sp 
Sp 



7a||7/3||n CT 
7a||73||75n (T 

7m||7HI7p||7<5||IT (T 
7u||7HI7p||7<5||75II, 



2A Q/3 , 
= 2<j(p a/3 , 

2(A^ u A p s + A^sA up — Afj, p A u s) 
= 2a(Ap V {p p6 + Apgip^u), 



(A.14 

(A.15 

(A.16 

(A.17 
(A.18 



In what follows, A 



{tptp)p U and A 



(<p<p)p U , where </V = F^/B and (p^ 



F pv /B are the tensor and dual tensor of the external electromagnetic field reduced to 
dimensionless form. The calculated nucleon spur can be reduced to 



Arcrcr 



2m p m n (A a/3 + icnp a(3 ) 



(g 2 a -l) + (l + g 2 a )(uAu) 



(A.19 
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ATCTCT 

JV a/3 



2m p m n 



;i - g 2 a )A aP + (1 + gfi [2{uA) a {uA)p - u 2 A a p 
-2g a a ((uA) a (0u)p + {(pu) a {uA)p 



(A.20 



Since tensors (jA.20 |) and (jA.19 |) have only longitudinal and transverse components, 
respectively, it is convenient to break up the lepton spur into the same structures. When 
it is convolved with the nucleon spur, only the completely longitudinal and completely 
transverse parts of the lepton spur will give a nonzero answer: 



L 



a/3 



M acT 
JV o/3 



-± 
J aP 



(A.21 



After simple transformations using the properties of the polarization operator (jA.12 J) , 
the corresponding constructions can be reduced to 



J a(3 



J af3 



2Sp 70||&||7 a ||p||(l -7 5 )IL 



2Sp 



7S-L&||7a±P||(l - 75)11. 



(A.22 
(A.23 



Using the above properties (|A.14 |) - ()A.18 |) . these expressions can be easily calculated and 
represented as 



L lp = ^[(P^)a(kA)^ + (kA) a (pA) p - A aP (pAk) 

+A({<pp) a {kA)p + {pk) a (Cpk)p), (A.24 ) 

= 4 (A a p - iif a p^j ((kAp) + . (A.25 ) 

For completeness, we present the convolution of the nucleon and lepton spurs with the 
origin in covariant form: 



ly aP 



J aP 



2(1+^) {{up)\\{uk){) + Ag a a{u{pp){u(pk) 



jvrcrcr 



L 



aP 



+2(1 + ag a ) 2 {(up)\\{uk)\\ + {up)\\{wpk) + {uk)\\{u(pp)) , (A.26 ) 
16m p m n (l + a) \(1 + g 2 a )(uAu) - (1 - g 2 a )] \(k~Ap) + (kipp)] . (A.27 ) 



In the selected frame of reference, where the nucleon medium is at rest [it = (1,0,0,0), 
(kp)\\ = (kAp) = k p — k 3 p 3 ], these expressions are 



Iy aP 



Arcrcr 
JV a/3 



J aP 



16m p m n (l + ag a ) 2 (p + p 3 )(k + k 3 
= 32m p m n (l + (r)g 2 a (p + p 3 )(k - h 



(A.28 ) 

J a p — ^"V" , «v- L ~r u jy a Kfo ~<~ fsjK^o — ">3J- (A. 29 ) 

Hence, it is easy to obtain the final expression for the square of the amplitude (jA.8 |) of 
the process being calculated: 

\M a \ 2 = 4(p + p 3 ) [(1 + °ga) 2 (k + k 3 ) + 2(1 + a)g 2 a (k - k 3 )] . (A.30 ) 

Note that in the case of nonrelativistic nucleons under consideration, the square of the 
amplitude depends only on the longitudinal electron and neutrino momentum compo- 
nents. 
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Appendix B 



The neutrino absorption coefficient 



Below, we give the details of our calculation of the absorption coefficient JC^ (JTHJ) for 
process (JTJ) involving neutrinos. We used standard expressions for an elements of phase 
volume: drS^ = {dp 2 dpf 1 L y L z ) / {2tt) 2 for charged particles and dn = (d 3 pV) / (2tt) 3 for 
uncharged particles. The proton distribution function was assumed to be the Boltzmann 
one. Eliminating the integration over the proton momentum due to the 5^ 3 ^-function, we 
note that the absorption coefficient does not depend on the P2 component of the electron 
momentum. Since P2 = eBx c in a magnetic field (x c is the centroid coordinate of the 
distribution of the electron wave function across the magnetic field) , we can eliminate the 
integral over the P2 component of the electron momentum: J dp2 = eBL x The remaining 
5-function in energy can be simplified. Since the electrons and neutrinos are ultrarela- 
tivistic particles and the nucleons are nonrelativistic particles, their momenta are of the 
order of P 2 ~ Tm (for nucleons) and p 2 ~ T 2 (for ultrarelativistic particles). Assuming 
the nucleon masses to be identical (where this does not cause any misunderstanding), we 
obtain 



(q) 2 (q\\ + k\\ - p\\) 2 eBn 



2m r 



2m r 



m. 



+ uj - e + (m n - m p 



A. 

2m r , 



eBn 



+ u - e + (ra„ - m p )J, 



fB.l 



where q± and gii are the particle momentum components across and along the magnetic 
field, respectively. Since the dimensionless parameter b = (eB) / (2m p T) < 1 we may 
neglect this term in one of the 5-function in the expression for the square of the S'-matrix 
element fll2jl by combining the two sums. In this case, the expression for the absorption 
coefficient is simplified and can be represented as 



"<2 „„ 2 



(2ny 



dy 



E 



l + exp[r] e -y]^n\ \b 



uj + m n -m p 
x o \ x — on — y -\ — 





fey /u) 


vn 


' X' 


I exp 


~b. 



f n F(ql,b)d s q, (B.2) 



Here, x = q\/(2m n T) and y = e/T. Note that the function F(q 2 ± ,b) has a simple 
asymptotic behavior for b — > 0, bn is a finite number. The latter condition follows from 
the fact that as the magnetic field weakens, the number of Landau levels that contribute 
to the sum under consideration increases in inverse proportion to field strength. In this 
case, the sum over the proton Landau levels may be replaced by an integral. It can also 
be shown that in this limit, 



lim 



1 

n\b 



cxp 



5(x — bn), 



(B.3) 
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where 5 is the Dirac delta function. Hence, it is easy to find that 

lim [F(ql, b)} = -. ( B. I 

1 + exp ( rj e — [to + m n — m p ) /T 



Since the function F(qj_,b) does no depend on q\ in the limit under consideration, the 
integration in (|B.2 |) can be brought to the end. Given the definition of the neutron 
number density, N n = 2/(27r) 3 (J f n d 3 q) we obtain the final expression for the neutrino 
absorption coefficient in the limit considered: 

KM = Gl™^ eBNn I'^Kl'-W" (B . 5 , 

~" 1 + exp (r] e - (u + m n - m p )/TJ 

The absorption coefficient for process © involving antineutrinos can be calculated in a 
similar fashion to give 

= ^^BN, d + (l-^W" (B.6 ) 

2n 1 + exp ( - r] e - (u - m n + m p )/T\ 

Appendix C The force density along the magnetic 

field for the neutrino scattering by nu- 
cleons 

In the low-energy limit, the effective Lagrangian for the neutrino scattering by nucleons 
is [5] 

C= ^(^(P0Ta(c„ + c a 75)M^))(^(^)T«(l + T5)^)). (C.l ) 

Here, Un and U u are the Dirac nucleon and neutrino bispinors; and c v and c a are the 
vector and axial constants of a neutral nucleon current. At low energies [22], c v = —1/2 
and c a ~ —0.91/2 for neutrons and c v = 0.07/2, c a ~ 1.09/2 for protons. Using the 
nonrelativistic density matrix for the polarized nucleons, p = m^(l +£<?), where ^ = 
±B/B is the nucleon polarization vector along the magnetic field and a are the Pauli 
matrices, it is easy to obtain the following expression for the square of the S'-matrix 
element for the neutrino scattering by nucleons per unit time: 



\S if \l _ (2tt) 4 G^ 



f 5^(V + k — V' — k') (cl + ScDuuj' + (cl - cl){kk') 



T 2V 3 luu> 

+2c v c a {S + S')(uk[ { + u'k l} ) - 2c 2 a {S - S')(uk\ - 

+ {c 2 v - cl)SS'(uu' + ikk')) + 4c^5"A;||A;[|l . (C.2 ) 

Here, k = (co,k) and k! = (u r , k') are the 4- momenta of the initial and final neutrinos; 
k\\ and hi are the components of these momenta along the magnetic field; and S, S' = 
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±1 (S = (£B)/B are the components of the polarization vectors for the initial and 
final nucleons along the magnetic field). The square of the S'-matrix element for the 
antineutrino scattering by nucleons can be derived from (jC.2 |) by the change k <-> k'\ 



\S„ 



if\i> 



\S if \l(k^k') 



(C.3 



Since we consider the neutrino scattering by nucleons in a strong magnetic field, for 
the subsequent calculations, it is necessary to take into account the contribution to the 
nucleon energy from the interaction of their magnetic moment with the magnetic field: 
Em = rriN + V 2 / '(2mjv) — (gNSeB)/ (2m n) where qn is the nucleon magnetic factor (g n ~ 
—1.91 for the neutron and g p ~ 2.79 for the proton). 

Next, we can simplify the 5-function in energy. Using the fact that V,V ~ y/m^T 
and k,k' ~ T and disregarding all terms ~ yjT /m^, we obtain 

/ , w - / , V 2 , a _„ eB (V + k + k^ 2 

5 (uj - J + E N - E N ) =6\u-u'+ g N (S - 5"); 



2m 



N 



2m N 



2m 



N 



8 I to — uo' — qn(S — S') 



eB 



2m 



N 



(C.4 



It is convenient to perform the subsequent calculations by separating out the contributions 
from S = S' and S = —S'. Eliminating the 5-function in momentum of the final nucleon 
and given that all terms linear in k\\ and fcf, give no contribution to the momentum transfer 
to the medium along the magnetic field, we derive the following expression for the density 
of the force acting on a unit volume of the medium per unit time: 



S=S' 



E 



S=S' 



S=-S' 



2Gl Co 
(2tt) 8 

s 



d S Vf N (V) / d s kf u (k) / d A k'(l - f v (k')) 



x-^—^u'kl - uk\ 2 )5 (u - u 1 ) , 

d 3 Vf N (V) j d?kf v (k) j d 3 k'(l-f u (k>)) 



UJUJ' 

2G 2 F c 2 a 
(2vr) 8 
S 

x (uj'kf\ + uk\ 2 )5 (uj - J - g N SeB/m N ) . 

UJUJ 1 11 



(C.5) 



(C.6 



Here, we use the fact that the nucleon gas is a Boltzmann one and, hence, 1 — /jv — 1. 
Since eB / (m^T) 1, the distribution function for the initial nucleons can be written 

as 



f N (V) ~ 1 + g N SeB/(2m N T) exp - V 2 /(2m N T) + Vn 



(C.7) 

where rj^ is the nonrelativistic nucleon degeneracy parameter. Interacting over the initial- 
nucleon momentum and using the fact that the nucleon number density is 



N 



N 



2/(2tt) 



exp 



V 2 /(2m N T) + VN 



d 3 V, 



(C.8) 
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we obtain the following expression for (|C5 J) and (jC.6 |) 

3f = ^^T 6 S f 1 + g N S-^-] h, (C.9 
11 s=s> (27r) 5 \ 2m N T J 

OO OO 1 2-7T 1 2tt 

h = J y 2 dy J y' 2 dy' J d X J d<p J d X ' J dtp' f u (y, X )(l - W,x')) 

-10-10 

x (V (1 - x 2 ) cos 2 ip - y(l - x' 2 ) cos 2 </) 5 (y - </) , (CIO 
G ^,(l + ^^)/ 2 , (C.11 



S=-S' (27r) 5 ' 

2tt 1 2?r 



J. 



2 







y 2 dy J y' 2 dy' J d X J dip J d X ' J d^f u (y, X )(l-My',x')) 

-10-10 

x (V (1 - x 2 ) cos 2 <p + y(l- x' 2 ) cos 2 <pi)s (y - y> - 9nS^\ . (C.12 ) 

Here, y = oj/T, y' = to'/T, x = cos 9, and x' — cos9', where 6 and 9' are, respectively, the 
angles between the momenta of the initial and final neutrinos and the radial direction; 
and ip and ip' are the azimuthal angles. Since the square of the S"-matrix element for the 
process involving antineutrinos can be derived by the change k <-> k', it is easy to note 
that the expression for the momentum transferred along the magnetic field in the process 
of antineutrino scattering by nucleons is given by the formal change c 2 — * — c 2 : 

^ = 9f(cS -►-<£). (C.13) 

For the integrals (jC.10 j) and (jC.12 I) to be calculated, we must use an explicit neutrino 
distribution function. We choose it in the form (JTTJ): f u (k) = <& v {r, x)F u (uj) where F u (u) = 
(1 + exp[u;/T — ?7y]) _1 , rj v is the neutrino degeneracy parameter. 

It is convenient to express the final result in terms of the mean neutrino-flux param- 
eters (J221), (|23|) . and (JHUj) . For the distribution function used, these parameters can be 
represented as 

J v = (47T)- 1 j * v (r,x)dn (C.14) 
(Xl) = (/ X 2 Mr,x)d^ (j ^(r, X )dn\ (C.15 ) 

(u?) = (J cu n+2 F v (uj)d 3 u?j (J u 2 F u (u)d 3 u^j (C.16 ) 

For completeness, we present the result of our calculation of the contributions to the 
force density from different nucleon polarizations separately: 

E 9fr ,|^ = _^^£_ NM ^ _ 1/3)i (c . 17 , 

s 
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-(<^> - 5T(ul)) (5/3 - (xl)) + 2J,(K 3 ) - hTM)) (l " 0£>) 



. (C.18 ) 



It is easy to verify that each of these expressions is zero for a thermal equilibrium of the 
neutrinos with the medium (T v — T, J u = 1, (xl) = 1/3), as follows from fundamental 
physical principles. 
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